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AN OPTIMAL APPROXIMATION OF ROSENBLATT SHEET BY 
MULTIPLE WIENER INTEGRALS* 
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Abstract. Let be the Rosenblatt sheet with the representation 

f [ [ f Q°‘{t,yi,y2)Q^is,ui,U2)B{dyi,dui)B{dy2,du2) 

Jo Jo Jo Jo 

where R is a Brownian sheet, \ < a, jJ < 1, (3“ and are the given kernel. In 
this paper, we contruct multiple Wiener integrals of the form 

f [ [ [ [ki{.yi,y 2 )~^°‘{ui,u 2 )~^^+ k 2 {yiVy 2 )^°‘{yi J\y 2 )~^°‘\yi-y 2 U~^ 

Jo Jo Jo Jo 

■ {ui V U2)^^{ui A - U2f~^]B{dyi, dui)B(dy2, du2), fci,^ > 0 , 

and obtain an optimal approximation of s). 


1. Introduction 

Self-similar processes are stochastic processes that are invariant in distribution 
under a suitable scaling of time and space. This property is crucial in applications 
such as network traffic analysis, mathematical hnance, astrophysics, hydrology and 
image processing. For this reason, their analysis has long constituted an important 
research direction in probability theory. The Hermite process is an interesting class 
of self-similar processes with long range dependence, it is given as limits of the so 
called Non-Central Limit Theorem studied in Dobrushin and Major j6], Taqqu HU. 
Let us briefly recall the general context. 

Denote by Hj{x) the Hermite polynomial of order j dehned by 

/ , , \ T 3d. 

Hjix) = i-iye 2 —e 2 , j = 1 , 2 ,... 

with Hq{x) = 1, and let the Borel function : M —)■ M satisfy E{g{^o)) = 0, < 

00 and 

00 ^ 

9{x) = Cj = -E[g{^o)Hj{fo)]- 

j=o 

The Hermite rank of g is dehned by 

k = min{j : Cj 7^ 0}, 

Clearly, k>l since E[g{^o)] = 0. 
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Let g he a. function of Hermite rank k and let (^ri)neN be a stationary centered 
Gaussian sequence with = 1 which exhibits long range dependence in the 

sense that the correlation function satishes 


( 1 . 1 ) 


r{n) := E{^o^n) =n »= L{n), 


where /c > 1 is an integer, H G (|, 1) and L is a slowly varying function at inhnity. 
Then, the Non Central Limit Theorem implies that the stochastic processes 


[nt] 


n 


H 




i=i 


converges, as n —>■ oo, in the sense of hnite dimensional distributions to the process 


( 1 , 2 ) 





Jo 




-(W- 


dsdB{yi)...dB{yk), 


\j=i 


where x+ = max{x,0} and the above integral is a Wiener-Ito multiple integral 
with respect to the standard Brownian motion {B{y))y^]s_ excluding the diagonals 
{Ui = yj}c ii k) is a positive normalization constant depending only on H 
and k such that E{Z^{1)Y = 1. The process {Z^{t))t>Q is called as the Hermite 
process of order k, it is H self-similar and has stationary increments. The class of 
Hermite processes includes fractional Brownian motion {k = 1) which is the only 
Gaussian process in this class. Their practical aspects are striking: they provide 
a wide class of processes from which to model long memory, selfsimilarity, and 
Holder-regularity, allowing significant deviation from fractional Brownian motion 
and other Gaussian processes. Since they are non-Gaussian and self-similar with 
stationary increments, the Hermite processes can also be an input in models where 
self-similarity is observed in empirical data which appears to be non-Gaussian. For 


k > 2, the process (1.2) is not Gaussian. When k = 2, the process (1.2) is known 


as the Rosenblatt process (see Taqqu [IE]). More works for the Hermite process 
and Rosenblatt process can be found in Albin [2], Leonenko and Aim |H|, Abry and 
Pipiras jl], Maejima and Tudor [9], Tudor in, Ghronopoulou et al. [5], Tudor and 
Viens pni. Torres and Tudor [H], Shieh and Xiao [15], Pipiras and Taqqu in, 
Maejima and Tudor [ini[II], Chen, Sun and Yan [1], Garzon, Torres and Tudor [7], 
Tudor [21], Yan, Li and Wu [22], Shen, Yin and Zhu [H] and the reference therein. 

Motivated by all these results, in this paper, we will prove the optimal approxi¬ 
mation theorem of Rosenblatt sheet based on the multiple Wiener integrals of form 


ffff 


I a—1 


[ki{yi,y2) ^“(mi,M 2 ) + A: 2 ( 2 /i V |/ 2 ) 2 “( 2 /i A 1 / 2 ) ^“| 2 /i-|/ 2 | 

■ {ui V U 2 )^y{ui A U 2 )~^y\ui - U 2 \y~^]B{dyi, dui)B{dy 2 , du 2 ), t,s > 0 
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with ki,k 2 > 0. Recall that Rosenblatt sheet with parameter ^ < a,/3 < 1 admits 
an integral representation of the form (see Tudor [21]), for t G [0,T], s G [0, S] 





t ps pt ps 



Qa(t, yi,y2)Qp{s, ui, U2)B{dyi, dui)B{dy2, du2) 


>0 JO Jo Jo 

where R is a standard Brownian sheet and K is the deterministic kernel given by 


( 1 . 3 ) 


with ch = 


= chS 2 ^ / {u — s)^ 2du for t > s, 


B{2-2H,H-J) 


R(-, •) represents the Beta function, H' = hM a^d dn = 


^ normalizing constant. Denote 

2/1, 2 / 2 ) = d//l[o,i](2/i)l[o,t](2/2 


' dK^' dK^' 

x,yi)^ — {x,y2)dx. 


'yi'Jy2 


dx 


dx 


In general, for every Borel measurable function ( G L^([0,T]^ x [0,5]^) the sto¬ 
chastic integral 


t ps pt ps 


MUi) ■■= 



((!/i,!/2,'‘‘i,U2)B(dpi,dui)B(d!/2,du2),t E | 0 ,T],s e | 0 ,S] 


'0 Jo Jo Jo 

is well-dehned, and the optimal approximation problem is to estimate 


(1.4) 


inf sup E(Z'^’^(t,s) — Mts(C))‘^- 

CeL2([o,T]2x[o,5]2)ig[o,r],ne[o,S] 


Noting that if the above minimum is attained at the function (*, then (C* > 0 a.e. 
In fact, we have 

£’[Z"’^(f,s) - Mt,s(C)]^ =2 / III C{yi2y2,ui,U2)du2dy2duidyi 


t ps pt ps 



Jo Jo Jo Jo 
yi,y2)Q^{s, Ml, U2)C{yi, 2/2, mi, U2)du2dy2duidyi 


0 Jo Jo Jo 


for all t, s > 0. 

If C( 2 /i, 2 / 2 ,Mi,M 2 ) < 0, then 

sup R[Z“’/'(f,s) -Mi,,(C)]2 > sup R[Z“-^(t,s) -Mi,,(|C|)]2 . 
ie[o,T],se[o,s] te[o,r],se[o,s] 

This gives the contradiction. Hence, we can assume that ki,k 2 > 0 and study the 
optimal approximation problem 


( 1 . 6 ) 


inf sup £i|Z‘‘''’((,s) - M,.,(C)]^ 

ie[o,T],se[o,5] 
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where 

^ = {C(2/1,I/2,Mi,M 2) =fcl(2/l,|/2)"5"(Mi,M2)"^^ + k2{yi V I/ 2 ) 5“(|/l A I/ 2 )"5"|- 2 / 2 !° 

■ (Mi V M2)2^(Mi a - M2|^“\ Aji, ^2 > 0}, 

since l/i, 2 / 2 ) < C'H,r{(|/i, 1 / 2 )“^^ + (l/i V 1 / 2 )^'^(|/i A |/ 2 )“^^||/i - 1/2|'^"H- For 
A; G /C, denote 

/(t, s, fci, k2) := 2E[Z“'^(t, s) - Mt,,(C)]^ t, s > 0. 

The similar approximation for the fractional Brownian motion and Rosenblatt pro¬ 
cess are hrst considered by Banna and Mishura [3], Mishura and Banna [12] and 
Yan, Li and Wu [22], respectively. 

The rest of this paper is organized as follows. Section 2 give the representation 
of the function /(t, s, fci, ^ 2 ) = 2E(Z°‘’^{t, s) — Mt^s{QY for C ^ In Section 3, we 
consider the function sup^gjg 7 .] 5 ] /(A, s, fci, ^ 2 ) in the compact rectangle interval 

|o,r] X [0,5']. In Section 4 and Section 5, we consider the optimal approximation 
in the two case A < 0 and A > 0, respectively. Two special cases be considered in 
Section 6. 


2. The representation of /(t, s, fci, ^2) 

In this section, we will give the representation of /(f, s, /ci, ^ 2 ) = 2E[Z^'^{t, s) — 
m,,,(C)]2 for c e /C. 


/f 7 2 

a{k 2 ) = 1 + ^B {1 -a,2a- 1)R(1 -(3,2(3-!)- ^C 2 {a)C 2 {( 3 ), 


Theorem 2.1. Let 

8 h 

a /3 ’ / V r- 5 r- ' a (3 

h{k 2 ) = C'i(a)C'i(/3) - 4^25(1 - a, a)B{l -/?,/?), 
where C'i(a) = daC^/B"^/! — a, |a), and 


1 rs 


C 2 (a) = dacl, 



r-“(l - s)5“-^(l - r)2“-^(s - r)“-Mrds, 


0 JO 


for all ki, k 2 > 0 and ^ < a, (3 < 1. Then we have 

^2-2a 2-2/3 

f{t,s,ki,k 2 ) = a{k 2 )t‘^°‘s‘^^ - 8 kib{k 2 )ts + 4kl- -r^--—, t G [0,T],s G [0,5']. 

(i — — (3)^ 

Remark 1. As an immediate result we have 0 (^ 2 ) > 0 and 

1 


( 2 , 1 ) 


b3h)(i - a)4i - < -Ah) 


for all a, (3 E (|,1) , since /(f, s, fci, ^ 2 ) > 0. Notice that 0 (^ 2 ) is also a quadratic 
equation in k 2 , we get 

^Cl{a)Cl{P) < B{1 -a,2a- l)B{l - f3, 2(3 - 1) 
for all \ < a, (3 < 1. 

2. Using the constant Ci, C 2 we give the main results and at the end of this paper 
we give the numerical of these constants (see Figure 1,2,3, f.) 
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Proof. It is easy to calculate that 

ot pS pt PS 



0 -'0 -'0 -'0 


t pt 



S PS 


0 ^0 

- d 


Q°'{t,yi,y 2 ){yi,y 2 ) '^°‘dyidy 2 

'*t pt pt 



Q°'{t,yi,y2)Q'^{s,ui,U2){yiy2) ^^°‘{uiU2) ^^^du2dy2duidyi 

(5^(<S, Ui,U 2 )(ui,U 2 )~^^dUidU 2 

( 1 / 1 I/ 2 )■“«“(« - 1 / 1 )- y 2 y^'^~^dudyidy 2 
{uiU 2 )~^u^{u — ui)^^~^{u — U 2 )^^~^ dudyidy 2 



'0 00 

la-1/ 


'0 ^0 ^y\'\/y 2 

r*s PS PS 


X 


dyCpi 



'0 Jo Jui\/U 2 

^2 j ^2 td 2 i 


= daci,dycj,B'\l - a, - P, ^/3)ts = Ci{a)Ci{P)ts, 


and 


t PS pt PS 



0 JO Jo Jo 


Q"(t,|/I,|/2)Q^(S,U1,M2)(|/1 Vl/2)^“(|/l Ays) -I/ 2 I 

■ {ui V M2)^^(mi a M2)“^^|mi - U2f~^du2dy2duidyi 


a —1 


t PU PU 


— ri 



0 Jo Jo 


X 


(I/1I/2) -1/1)2“ ^ 

■ (yi V |/2)5“(|/i a |/ 2 )"^“||/i - y 2 \‘^~^dyidy 2 du\ 

r riuiu 2 )--^^u^{u-u,)-^^-\u-u 2 )-^^-^ 

Jo Jo Jo 

■ (til V U2)^^{ui A U2)~'^^\ui — U2\^~^duidu2du] 


^ C 2 {o) ^ 2 a ^ 2 {P) ^ 2 y 
a P 

Then for all t G [0,T], s G [0, S], we have, 

nt PS pt PS 



0 JO JO JO 

On the other hand, 

pt PS pt PS 


yi, y2)Q'^{s, ui, U2)C{yi, y2, Ui,u2)du2dy2duidyi 

hCPa)CPP)tS + fc2 ^^(°')^2(/^) /2a^2/3_ 

ap 



0 -'0 JO JO 


C {yi,y2, Ui, U2)du2dy2duidyi 

t ps pt ps 



^ 1 ( 1 / 11 / 2 ) “(M 1 M 2 ) ^du 2 dy 2 duidy 2 


'0 JO JO JO 

pt PS pt PS 


+ 



^2(1/1 V 1/2)“(1/1 A 1/2) "I1/1-I/2P" ^ 

'0 JO JO JO 

{ui V U 2 YPJ 1 A n2)“^|Mi - U2\‘^'^~‘^du2dy2duidy2 
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s s 

+ 2 /CiA; 2 (|/i|/ 2 )" 5 “(|/i V|/ 2 ) 5 “(|/i A|/ 2 )" 5 “||/i-1/2I 

Jo Jo Jo Jo 


a—l 


■ {U 1 U 2 ) 2^(mi V M2)2^(mi A M 2 ) 2/^|mi-M2|^ ^du2dy2duidy2 

'■= I 1 + I 2 + h, 


where 


t pt 


h = kl / {ym) ‘^dyidy 2 / / (M 1 M 2 ) ^duidu 2 

Jo Jo Jo Jo 


^2-2a^2-2/3 

(l-ani-fjy 


kl 



h = kl{[ [ (yi V|/2)“(|/i A 1 / 2 ) “li/i-|/2r“ ^dyidy2) 

Jo Jo 

f*s ps 

(mi V M2)^(|/i a M2)“^|mi - U2\‘^'^~‘^duidui) 

'0 JO 

r p ry2 rt ft 

= kl y 2 yi''ijj 2 -yif'"~^dyidy 2 + j / yiy 2 ’^{yi-y 2 f'^~^dyidy 2 \ 

_J 0^0 J 0 J y2 

s pU 2 pu ns 

/ U 2 U{^ {U 2 — duidu 2 + / / U^U 2 ^ {ui — U 2 Y^~‘^ duidu 2 

IJ 0 Jo Jo Ju2 

= kl{2 [ yl'^~^dy2 j — uY'^~‘^du) ■ {2 f u\^~^du2 f u~^{l — u)‘^^~'^du) 

Jo Jo Jo Jo 


= \b{ 1 -a,2a- 1)5(1 - l3,2/3 - 

ap 


h = 2kik2{ [ [ ( 1 / 1 I/ 2 ) ^“( 1/1 V |/ 2 )^“(l/i A 1 / 2 ) - 1/21“ ^dyidy 2 ) 

Jo Jo 

ns ns 

■(/ / {UiU 2 )~^^{Ui\/U 2 )^'^{Ui AU 2 )~^^\Ui - U 2 f~^dUidU 2 ) 

Jo Jo 

ft fy2 ft ft 

= 2kik2{ / yi'^iy2-yiT~^dyidy2+ / y^'^iyi - y2T~^dyidy2) 
J 0 J 0 J 0 J y2 

ns PU2 PS ns 

■ ( / / m)~^(m2 — Ui)^~^duidu 2 + / U 2 ^{ui — U 2 )^~^dUidU 2 ) 

J 0 Jo J 0 J U2 

nt nl ns nl 

= 2 kik 2{2 / dyi / u~°'\l — u\°'~^du) ■ {2 / dui / u~^\l — u\^~^du) 

Jo Jo Jo Jo 

= 8/ci/c25(1 — a,a)B{l — (3, f3)ts. 


Hence, 

pt ps pt ps 



dyi,y 2 ,ui,U 2 )du 2 dy 2 duidyi = 


kl 


0 Jo Jo Jo 


(1 - a)2(l -/3)' 


^2-2a^2-2l3 


kn 


+ ^5(1 - a, 2a - 1)5(1 - /3, 2/3 - 1)^^^^ + 8A;iA;25(l - a, a)5(l - /3, f3)ts, 
ap 
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for all C E 1C. It follows that 

fit, s, h, k2) = 5) - Mt,siC)f 

Ciyi, 1/2, Ui,U2)du2dy2duidyi 
(5“(t, yi, y2)Q'^is, Ml, M2 )C(//i, 1/2, U2)du2dy2duidyi 

^2-2a^2-2l3 


/*S pt PS 

^ ^2q:^ 2/3 _j_ ^ I I I I /-2, 



>0 Jo Jo Jo 

pt PS pt PS 



0 ^0 ^0 ^0 


= a(/c2)t^"s^^ — 8kib{k2)ts + 4kl 


(l-a)2(l-/3)2' 


This completes the proof. 


□ 


3. The maximum value of fit,s,ki,k2) 

In this section, in order to obtain the optimal approximation with ki,k2 > 0. 
We need to hnd the maximum value point Po(^o,'So) in Ihe open rectangle interval 
(0,T) X (0, S), and the maximum value point at the boundary of rectangle interval 
|o,r] X [0,5']. Thus, we can get the sup fit,s,ki,k 2 ). 

t&[0,T],s&[0,S] 

Lemma 3.1. The function fit, s, ki.k 2 ) at the open rectangle interval (0, T) x (0, S) 
can’t get the maximum value point Po(^o,'So), to G (0,T),So G (0,5). 


Proof. If the function fit, s, ki, k 2 ) have maximum value point Po(to, Sq), to G (0, T),sq G 
(0,5), then Po(to,So) must be the stagnation point of the function 

(t,s) HA fit,s,ki,k2). 


That is. 


dt 

Solving equations set: 


(9/(t,s,/ci,/c2)| dfit,s,ki,k2)^ 

l-Po(io,5o) *^7 


dfit,sM,k 2 ) ^ 2tl-2«52-2/3^^^^^^)^4a-2^4/3-2 

- 4kibik2)t^‘^-h^^-^ + 


4kl 


os 

- 4kibik2)f‘"-^s^^-^ + 


(l-a)(l-/3)2 

Akl 


= 0 , 


(1 - a)2(l -13)' 
Let X = elementary calculation can obtain 

ki 


= 0. 


X = 


and 


6(A;2)(l-a)2(l-/3)2’ 
'/af^ = 26(A;2)(1 - q;)(1 - /3). 
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Hence, we have 

ki 

/(to, So, fci, ^ 2 ) = toSo[a(/c 2 ) ,... -^ - 4/ci6(fc2)] = 0. 

b{k2)[l - a)^{l - fjy 

But /(£,£, 0,0) = > 0, for any £ > 0. That is contradiction. Therefore, the func¬ 

tion fit, s, ki, ^ 2 ) at the open rectangle interval(0, T) x (0, S) can’t get the maximum 
value point Po(^o,So), h G (0,T),so ^ (0,5'). □ 


At the boundary 0 x [0,5'], and [0,T] x 0, we have /(t, s, fci,/C 2 ) = 0. So, we 
need to consider the case of extremum value point on the boundary T x [0, 5'] and 
[0, T] X S, with ki,k 2 > 0. We only think about the case of the boundary T x |0,S]. 
In the same way, we can get case of the boundary [0, T] x S. 

It follows from Theorem 2.1, we know that 


rp2-2a 2-2/3 

(3.1) /(T, s, fci, k 2 ) = aik 2 )T^^s^^ - ^k^hik 2 )Ts + 

for all s G [0,5']. Differentiating (3.1) with respect to s leads to 

Skf 


(3.2) /,(T, s, ki, k 2 ) = 2 ( 3 aik 2 )T^‘^s^^-^ - 8 kibik 2 )T + 


(3.3) 


(l_tt)2(i-/3) 

Let /s(T, s, ki, ^ 2 ) = 0 and x = kis^~‘^^, which implies that 

42-2-2a^2 


2,2-20^1-2/3 


F{x) := /3a(fc2)T"" - Abik 2 )Tx + 


il-ani-f3) 


= 0 


and the discriminant A of the quadratic function Fix) is 


If A < 0, then 


A = 16T\ibik2))^ 


/3aik2) . 

il-a)yi-/3)^ 


sup fiT,s,ki,k 2 ) = fiT,S,ki,k 2 ) 

se[o,S] 

(3.4) rp2-2a 02-2/3 

= T^^S^f^aik2) - 8kMk2)TS + 4fc?--—-— 

(1 — a)^(l — p)^ 

Using the same method, on the boundary [0,T] x 5, we have 


sup fit, S, ki, k 2 ) = /(T, S, ki, k 2 ). 

i6[0,T] 


If A > 0, then the equation (3.3) has two real roots as follows 


Xi = 


il-a)\l-y) 

22"i-2q 


( 6 (^ 2 ) + 



I3aik2) . 
(l-a)2(l-/3)^’ 


and 


X2 


il-a)\l-y) 

2J^1-2q 


(^^2) 



/3a(/c2) . 

(l-a)2(l-/3)^’ 


1 1 1 1 

which says Si = kl°‘~^x^ 2 a-i ^ _ ^ 20 ,- 1 ^^ 2 a-i stagnation points of 

the function s 1 —>■ /(T, s, ki, /c 2 ). Hence, Si, S 2 are the points of local maximum and 
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minimum, respectively, since the monotonicity of the function s ha f{T,s,ki,k 2 ). 
This implies that 

sup /(T, s, fci, A;2) =/(T, S', fci,/c2), if si > S', 

se[o,s] 

sup /(T, s, fci, A; 2 ) = max{/(T, si,A;i,/c 2 ),/(T, S', A;i,/c 2 )}, if si < S. 
se[o,s] 

Using the same method, on the boundary [0,T] x S', we can find a ti, such that 


(3.5) 


(3.6) 


snp f{t,S,ki,k 2 ) = f{T,S,ki,k 2 ), if ti > T, 

ie[o,T] 

sup /(t,S,/ci,/c 2 ) = max{/(ti,S,/ci, A;2),/(7', S',/ci, A; 2 )}, if U < T. 

te[o,T] 


4. The optimal approximation, case A < 0 
Theorem 4.1. If A < 0, then we have 

'j-^2—2a c2—2/3 

inf sup f(t,s,k„h,) = r^“S^^a(k;)-Sk;b{k-)TS + 4{k;)\ - 

Ce/Crx[0,s] 

where 


I Q —1 


^ = {C( 2 /i,I/ 2 ,Mi,M 2 ) =fcj!'( 2 /i,l/ 2 ) 2 “(mi,m 2 ) + A:;(i/i V 3/2 )^“(|/i A 3/2) ^“| 2 /i-|/ 2 | 

■ (mi V U2)^^{ui A M2)~^^|mi - M2|^~\ ki, k2 > 0} 
and {kl, kf) is the stagnation point of the function 

{ki,k2) HA /(T, S', ki,k2), 

here 

(4.1) 

_4B{1 - a, 0)5(1 - /3, /3)C2{a)C2{/3) - 5(1 - o, 2o - 1)5(1 - /3,2/3 - l)C'i(o)C'i(/3) 
^ “ 8o/3(i - o)2(i - yysyi - o, 0)52(1 -y,y)- 5(1 -0,20-1)5(1 - /3,2/3 -1) 

X (1 - o) 2 (l - 

(4.2) 


k* = 


C'2 (o)U2 (/3) - 4o/3(l - o)2(l - /3)25(1 - o, o)5(l - /3, /3)C'i(o)C'i(/3) 

^ “ 5(1 - o, 2o - 1)5(1 - /3, 2/3 - 1) - 8o/3(l - o)2(l - yys^l - o, o)52(l - /3, /3)' 

Proof. It follows from equation (3.4), when A < 0, 

rp2-2a 02-20 

sup /(T, s, fci, k 2 ) = T^^S^^a{k 2 ) - ^hh{k 2 )TS + ^{kif- --— 

se[ 0 , 5 ] [l-ay[l-py 

for all ki,k 2 > 0. Let now {ky kf) be the stagnation point of the function 

{ki,k 2 ) HA /(T, S', fci,fc2)- 


An elementary calculation can obtain k\, which can be denoted by (4.1) and (4.2), 
and the Hessian matrix H as follows 


H 


d‘^f{T,SMk 2 ) d^f(T,S,kik 2 ) 
dkf dk\dk2 

d^f{T,S,kik 2 ) d^f(T,SMk 2 ) 
dk 2 dk\ dk^ 
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and 


d^f(T,S,kik2) 8T2-2“^2-2/3 


dkf (l-a)2(l-/3)2 

d^f{T,S,kik2) d^f{T,S,k,k2) 


> 0 , 


dkidk2 dk2dki 

= 32TSB{1 - a, a)B{l -/3,/3), 

d^f{T,S,kik2) 


dki 


a/3 


-5(1 -a,2a- 1)5(1 - (3,2/3-1). 


So, 

|H| = 64T"5 


2/3 1 ) - 1 ^ 02/1 ^ ^.\u2/ 


165^(1 — a, a)B^{l — (3, (3)) > 0 


a(l - a) 2 / 3 (l -/ 3)2 

for all ^ < a, /3 < 1, since 5^(1 — a,a) < ■ This means that the minimal 

valne of {ki,k 2 ) e-)■ f(T,S,ki,k 2 ) is achieved at point (/cjj',/;^). This completes the 
proof. □ 


5. The optimal approximation, case A > 0 


Lemma 5.1. //A > 0, then we have 

S,ikl,k* 2 )<S, 

and 

S 2 iki,k;)<s, 

where {k\, k^) is the stagnation point of the function {ki, k 2 ) e-)■ snp^gjQ^] /(T, s, ki, /C 2 ). 


Proof. We split the proof in two steps. 
Step one. It is easy to obtain that. 


Ti + a;2 _ b{k 2 ){l - a;)^(l - /3) 

~ Ti-2«A;i ' 


bi bi *2 1 

since S 2 > Si. Hence, when ki = kl and ^2 = ^ 2 , 

(5.2) 

b{k* 2 ){l - afil - (3) {Ci{a)Ci{(3) - 4A;*5(1 - a, a)5(l - (3, /3))(1 - a)^{l - /3) 


Ti-2“A;^ 


Ti-2«A;* 


Si- 2 g 2 

> 


1 - /3 ^2/3-1 

for (3 G (|, !)• This proves that Si(/c^, k^) < S. 
Step two. 

df{T,s,klk *2 


ds 


= 2/3a{k;)T^^s^f^-^ - 8klb{k;)T + 


8 {kl)^ 


_rj^2-2a ^l-2l3 


(1 -a)2(l -/3)' 

= 2da(k*)T^'"s^^~^ - 3 ^(^2)(1 ~ ~ )^ y2-2a 1-2/3 

2 / 3 a{k^)I 8 yi_2,^i_2;3 b{k^)I + 

= 2/3T^'^s‘^^-\a{k;) - 4b‘^{k;){l - af{l - (3f), 
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Since, 


Hence, 


dfiT,s,kl,k*) 


a{k *2 




ds 


=s > 0. 


U=5 


Thus, S < Si or S > 32 - By step one, we know S' > Si, so 52 (^ 1 , < S. 

Lemma 5.2. Denote 


□ 


h{ki, fca) = f{T, Si, ki, fca) - f{T, S', fci, /ca)- 


Then the equation h{ki, ^ 2 ) = 0 has two solutions k[ and ki, which satisfy 0 < < 


ki, ^\k.=k: >0 =0- 


dh 


dki \ki=k-^ 


dki \ki=ki 


Proof. By Theorem 2.1 , we have 


f{T, si, ki, /C 2 ) = T^"s4a(/c2) - 8 kib{k 2 )Tsi + Ak'f 


T 


2 -2 q ^2-2/3 




n^ 2 - 2 a(ki\ 2 ^ 
.2 ^ 


2 - 2/3 


n 1 23 1 

= r-( 

2/3 

:= k^~^(p{k2) = kl^~^ip{k2). 

Let /i(fci, fca) = 0, which implies that 
(5.3) 

2/3 T-'2-2a C2-2/3 

kl^-^^{k 2 )=T^^S^^a{k 2 )- 8 kih{k 2 )TS + Akl- --— = /(T,^, hM)- 

(1 — a)^[l — /j)^ 

j^2-2a ^2-2/3 


Differentiating (|5.3|) with respect to ki 

1 


(5.4) 


23 


23-1 


k^^-^ip{k2) = -Sb{k2)TS + Ski 


(l-a)2(l-/3)2' 


multiplying by ( 2 ^-^)fci both side of (5.3) leads to 

kp-Mh) = 


rp2-2a ^2-2/3 


(2/5 - l)ki 


2/5-1 


2 / 1-1 
It follows that 

rp2-2a 02-2/3 

23 a{k 2 )T^^S^^ lQ3b{k2)TS 


(2/5-l)(l-«)2(l-/5)2- 


(2/1 - l)fci 


2/5-1 


+ 8/lfci 


j^2-2a ^2-2/3 


(2/l-l)(l-a)2(l-/l)2- 


This implies that 
(5.6) 

(4/3-4)r=-^“S^-^'’fcJ+45(*:2)rS(l-a)“(l-/3)^*:,-/3a(*:2)(l-a)=(l-/3)T^“S^'’ =0. 
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This is a quadratic equation in ki with the two roots 

^ 2(1 - /?)Ti-2«5i-2/3 

fe(fc2)(l-a)2(l-/3)2- 
— 2(1 - /?)Ti-2«5i-2^ 

since 

D = b‘^{k 2 ){l — a)‘^(l — — /5a(/c2)(l — /9)^(1 — > 0. 

It is easy to check that ki is the solution to the equation set 

f h{ki, k 2 ) = 0, 

I §-Akuh) = o. 

In the follows, we will prove that ki is not the solution of the equation h{ki, ^ 2 ) = 0. 
In fact, 

h{h, k 2 ) = f{T, Si, h, k 2 ) - f{T, S', h, k 2 ) 

2/3 Arp2 — 2a 

= xl^-^^{k2)S^^ - a{k2)T^^S^^ + 8b{k2)TS^^X2 - 

2 P Arp2—2a 

= - a{k,)T-‘‘ + S,h{k2)Tx2 - (1 _„'2n 

= 'T^S^‘‘(C-)^'{iKk 2 )Txi - a(k,)r-‘) - (ib(k^)Tx^ - a(k,)T^‘‘)), 

1 — p Xi 


since 


2/3 

ip{k2) = a{k2)T^"x^ - 8b{k2)xf-^ + 


4T 


2-2a 


(l-a)2(l-/?) 


2^1 


2/3-2 


By (3.3), we have 


Xi + X2 = 6(fc2)T^" ^(1 - a)^(l - /d), 2:1X2 = -/3a(A:2)T'^“ ^(1 - a)^(l - /d), 

since xi,X 2 are the root of equation f]a{k 2 )T'^°‘ — 4 h{k 2 )Tx + = 0. Thus, 

Ab{k2)Txi _ X 1 + X 2 Ab{k2)Tx2 _ + X 2 

~ P I > \rn2a ~ P~ 


a(A:2)T2« 

Let X = ^ e (0,1). Then 


X 2 


a{k 2 )T 


Xi 


;^)^(46(A:2)Txi - a(A:2)T2“) - (46(fc2)Tx2 - a{k2)T‘^^) 

Xi 


2o,.X2.^./3(Xi +X 2 ) .. /^(Xi+Xa) 


= a(k 2 )T^^((—)^( 

Xi X2 

2/3 


-1)- 


Xi 


- 1 )) 


= a(k 2 )T^‘"(x^ (/3x-^ + /3-l)-/3-/3x + l) 

:= a(k2)T^'^p(x). 
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It is easy to check 0(x) > 0. In fact, 0(0) = 1-/3, 0(1) = 

0'(1) = O 


0"(x) 


20(1-0) 

(20 - 1)2 


2-20 1-20 

a: 2/3-1 (a: 2/3-1 - 1) > 0 


0, and 0'(O) 


- 0 , 


for all X G (0,1) since 20 > 1. This shows that the fnnction 0(a:) is convex on (0,1) 
and 4> {x) is increasing strictly on (0,1), which implies 4> {x) < 0. It follows that 
0'(x) is decreasing strictly on (0,1) and 


0 (a;) > 0(1) = 0 


for all X G (0,1), Thus, h{ki, ^ 2 ) > 0. 

On the other hand, h(0, ^ 2 ) = —a{k 2 )T‘^°‘< 0, it follows that the equation 

h{ki, ^ 2 ) = f{T, si, fci, ^ 2 ) - f{T, S, ki, /C 2 ) = 0 

admits a root, denoted by k[, on (0,/ci). Noting that the ki ha /(T, si,/ci, ^ 2 ) is 
convex and increasing, we hnd that the equation 


h{ki, ^ 2 ) = f(T, si, ki, ^ 2 ) - f(T, S, ki, /C 2 ) = 0 

admits two roots at most since the function ki ha f{T,S,ki,k 2 ) is a quadratic 
function. Thus, k[ is unique in (0, fci) and >0. □ 

Now, we consider sup^g[o,s] /(^) ^ 2 ) at the case A > 0. 

Theorem 5.1. //A > 0, k[ is given in Lemma 5.2. We have 

max{/(T, si, fci, ^2), /(T, S', fci, /C2)} = /(T, S', fei, ^2). 


Proof. By Lemma 5.1, we have 


sup f{T,s,ki,k 2 ) 

s£[0,S] 

It follows from Lemma 15.21 that 


max{/(T, si, ki, /C 2 ), f{T, S', ki, fca)}- 


max{/(T, si, fci, ^ 2 ), /(T, S, ki, ^ 2 )} = /(T, si, fci, k2)l^i,,>k[} + fi'^y ^y k2){k,<k[}y 


which implies that 


max{/(T, si, ki, ^ 2 ), /(T, S', ki, fca)} = f{T, S', k[, ^ 2 ), 


since ki ha f{T,si,ki,k2) is increasing and /(T, S,/ci, ^2) is decreasing for ki < 
k[. □ 


Remark 3. By Theorem |5.1 and (3.5[), we have 


sup f{T,s,ki,k 2 ) = f{T,S,k[,k 2 ), 
se[o,5] 


where k[ G (0, fci). 

Use the same method, on the boundary [0,T] x S', there also exists a k'l > 0, such 
that 

sup f{t, S', ki, k 2 ) = f{T, S', kl, k 2 ). 
te[o,r] 
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From all above, we have 

sup f{t, s, ki, k 2 ) = max{/(T, S', /c^, ^ 2 ), /(T, S', k:[, ^ 2 )} 
te[o,r],se[o,S] 

Theorem 5.2. l)If f{T, S,k[,k 2 ) > f{T, S', /C 2 ). Then minimal value 

inf sup f{t,s,ki,k 2 ) 
ie[o,T],se[o,S] 

is achieved at point (T, S', k[, k' 2 ) and equals to f{T, S', A;^, A;^)- 
/(T, S', k'l, k 2 ) < f{T, S', k^, k 2 )- Then minimal value 


inf sup f(t,s,ki,k 2 ) 

. y-~ ± V \ J J ) / 

ie[o,T],se[o,s] 

is achieved at point (T, S', k'[^ k'^) and equals to f{T, S', k'l, k'l). 
Proof. 1). If /(T, S', k'^, k 2 ) > f{T, S', k'l, k 2 ), then 


sup f{t,s,ki,k 2 ) = f{T,S,k^,k 2 ) 
ie[o,T],se[o,s] 

= T‘^^S^^a{k2) - S}^^h{k2)TS + 4(A;' 


j^2-2a g2-2l3 

(1- 0)2(1 -/3)2 


42-2a^2/3 

' a(3 


5(1 -0,20- 1)5(1 -13,213- l))kl 


+ (32A;;T^5(1 - o, o)5(l -(3,13)- 

a(3 

+ - 8A;;T^C'i(o)C'i(/ 3) + _ ^j2(f_ 


This is a quadratic equation of A; 2 , and '^'^^J^^^ 5(l — o, 2o — 1)5(1 — (3,2(3 — 1) > 0. 
It is easy to find that, when 


^2 = 


C2{a)C2{(3) - Aa(3k'fr^-^^S^-^^B{(i - a, o)5(l - (3, (3) 
5(1 - o,2o - 1)5(1 - (3,2(3- 1) ^ 


we have 


inf sup /(t,s, A;i, A; 2 ) = inf/(T,^, A;i, A; 2 ) =/(T,^, A;i, A; 2 ). 
ie[o,T]se[o,S] 

2). If /(T, S', A;), k 2 ) < f{T, S', kl, k 2 ), we have 

inf sup f{t, s, kl, k 2 ) = inf/(T, S, kl, k 2 ) 
C6^te[o,r]se[o,s] CeK. 

= f{T,s,ki,k';), 


where 


C2{a)C2{(3) - Aa(3k'lT^-^^S^-^f^B{l - a, o)5(l - /3, /3) 
5(1 - o,2o - 1)5(1 - (3,2(3- 1) 


This completes the proof. 


□ 
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6. Two SPECIAL CASES 

In this section we consider two special classes of the approximation function ( E 1C 
First, if /c 2 = 0, we only consider the boundary T x |0,S]. 

Theorem 6.1. Let ICi = {C( 9 i>» 2 , > 0}. 

(1) IfCl(a)Ct(li) - < 0, then 

itif sup 2 £:(Z"’'’((,s)-A 4.(C))2 = Y“S2'’(l-4Cf(a)Cf(^)(l-a)^(l-.8)^) 

te[o,T],se[o,S] 

k - CL{a)CLmi-»ni-0r 

LULilL rv — j'l-2agl-20 

(2) IfCl{a)Cm - (i-g)?i -„)2 > 0, then 

inf sup 2 E{Z-'^{t, s) - MUOf = f(T, S, k\, 0) 

CeA2 t6[o,r],se[o,5] 

where C(i/i, ?/ 2 , mi, M 2 ) = A;)(i/i?/ 2 )“^“(miM 2 )~^^ and k'l is the smallest root of the 
equation f{T, si, k, 0) — /(T, S, k, 0) = 0. 

Proof. For ( E ICi we have 

2 E{Z^’f^{t,s) - MtAOf = f{t,s,k,0) 

= _ sk,C,(a)CMts + 

and Ai = 64A:^T^[Cf (a)Cf (/9) — which competes the proof. □ 

Second, we consider the case ki = 0. 

Theorem 6.2. Let 

^2 = {C(2 /i,I/2,Mi,'W2 ) = k{yi V y 2 )^°‘{yi Ai/2)"^“|2/i -i/ 2 |“"^ 

■ (mi V M2)^^(mi a M2)“^^|mi - U2\^~^,k > 0} 


we have 


with k* 


inf sup 2E{Z^’^{t,s) - Mt^si.0? 

C 6 Y 2 te[o,T],se[o, 5 ] 

^ _ ^Cl{a)Cl{ 0 ) _ , 

^ a 0 B{l-a, 2 a-l)B{l- 0,20 -ly 


c 2 {ot)C 2 (y) 

S(l-o,2o-l)B(l-/3,2/3-1) • 


Proof. By Theorem 2.1 


we have 


inf sup 2E(Z“’^(f,s)-Mi,,(C))^ 
C6Y2 tg[o,r],se[o,S] 


inf a(k)T'^^S^y 
C 6 A 2 


The function 


Ak'^ 

a{k) = l + —B{l 
a0 


a, 2a - l)B{l-0,20 


CL 

l)--Ce(a)Cm 
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is a quadratic equation in k, then 


inf a(A;)T2"^2/3 _ a(k*)T^^S^^ 

Ce^Ca 


= (1 


a/3B{l 


4C|(a)C|(/3) _ 

a, 2a - 1)B{1- 13,2/3 - 1)^ 


with k* 
and 


C2(«)C2(/1) 

_B(l-o,2o-l)B(l-/3,2/3-1) ‘ 


C(|/1,I/2,Mi,M2) 


C2{a)C2{(3){yi V y 2 )^°‘{yi A 1 / 2 ) - 1/2 

B{1-a,2a-l)B{l- 13,213-1) 

■ {ui V M2)^^(mi a n2)"^^|Ml - U2f~^ 


where 1 / 1 , 1/2 > 0 and ni,n 2 > 0. This completes the proof. □ 
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Figure 1. a ha Ci(a) 
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Figure 2. { a ,( 3 ) ha Ci { a ) Ci {( 3 ) 
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Figure 3. a h-)- C'2(a) 



Figure 4. { a ,( 3 ) ^ C2(a)C2(/3) 













